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Electron-cyclotron maser emission from white-dwarf pairs 
and white-dwarf planetary systems 
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ABSTRACT 

By analogy to Jovian radio emissions powered by the electromagnetic interaction 
between Jupiter and its moons, we propose that close magnetic-nonmagnetic white- 
dwarf pairs and white-dwarf planetary systems are strong radio sources. A simple 
model is developed to predict the flux densities of radio emission generated by a loss- 
cone-driven electron-cyclotron maser. The radio emission from these systems has high 
brightness temperatures, is highly polarized, and varies on a periodic cycle follow- 
ing the orbital rotation. Masers from magnetic-nonmagnetic white-dwarf pairs, with 
orbital periods < 10 min, are expected to be detectable over a wide range of radio fre- 
quencies. Terrestrial planets in close orbits about magnetic white dwarfs, with orbital 
periods < 30 hr, can also produce detectable radio emission, thus providing a means 
to identify Earth-sized extrasolar planets. 
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1 INTRODUCTION 

The operating requirements for electron-cyclotron masers 
are (i) a population inversion in the electron distribution 
and (ii) a magnetized plasma in which the electron-cyclotron 
frequency Q e exceeds the plasma frequency ll> p (e.g. Dulk 
1985). The first condition can be achieved when the mag- 
netic field geometry in the source region allows the devel- 
opment of an anisotropy in the electron distribution. One 
common example is the loss-cone electron distribution (e.g. 
Wu & Lee 1979), where an electron pitch-angle anisotropy 
develops within a magnetic flux tube with converging field 
lines at each footpoint. Large pitch angle electrons are mag- 
netically reflected, whereas small-pitch-angle electrons are 
lost through collisions with high density plasma at the foot 
of the flux tube. The second condition is satisfied in magne- 
tized plasmas with a relatively low electron density and/or 
high magnetic field strength. 

These conditions are satisfied in a variety of astronom- 



* E-mail: willes@physics.usyd.edu.au 
f E-mail: kw@mssl.ucl.ac.uk 



ical systems, including the solar corona (Melrose & Dulk 
1982; Robinson 1991), terrestrial and planetary magneto- 
spheres (Zarka 1992; Farrell, Desch & Zarka 1999: Ergun et 
al. 2000), and solar- like stars and magnetically-active binary 
stars (e.g. Melrose & Dulk 1982; Abada-Simon et al. 1994; 
Slee, Haynes & Wright 1994; Trigilio et al. 1998). Electron- 
cyclotron maser emission is characterized by high degrees 
of circular polarization (typically nearly 100%), high bright- 
ness temperatures (T ^> 10 s K), and narrow beaming. 

Planetary radio emissions are typically solar-wind 
driven, with the exception of a significant proportion of 
Jovian radio emissions, which are correlated with the or- 
bital phase of Io (Bigg 1964), and to a lesser extent with 
the phases of Ganymede and Callisto (Menietti et al. 1998, 
2001). In the unipolar- inductor model for the Jupiter-Io in- 
teraction (Piddington & Drake 1968; Goldreich & Lynden- 
Bell 1969) a current circuit is set up along the magnetic field 
lines connecting Io to Jupiter, closing in the Jovian iono- 
sphere. The circuit is driven by the e.m.f. induced across 
Io, which is a conducting body, as it traverses the Jovian 
magnetic field. The existence of a Jupiter-Io current circuit 
was confirmed by the detection at ultraviolet and infrared 
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wavelengths of a footprint in Jupiter's upper atmosphere, 
which remains fixed at the position of the footpoint of the 
Io flux tube as Jupiter rotates (Connerney et al. 1993; Clarke 
et al. f996). The unipolar- inductor model predicts the ex- 
istence of point-like auroral footprints. The observed tail 
of ultraviolet emission, extending up to 100° in longitude 
downstream of the bright Io footprint (Clarke et al. 2002), 
reveals that plasma-inertia effects are also important in the 
Jupiter-Io interaction, in addition to the unipolar inductor. 
The plasma-inertia effects are due to the presence of a dense 
plasma torus about Io's orbit, consisting of ionized volcanic 
material from Io, which is "picked up" into corotation with 
the rotating Jovian magnetic field (Brown 1994). The au- 
roral ultraviolet footprints for the Jovian moons Ganymede 
and Europa appear as point-like sources (Clarke et al. 2002), 
consistent with the unipolar-inductor model, despite the 
fact that the Ganymede- Jupiter interaction is complicated 
by Ganymede's embedded magnetosphere (Kivelson et al. 
1996). An auroral footprint for Callisto has not been de- 
tected because the footpoint of the Callisto flux tube is co- 
incident with the main auroral oval (Clarke et al. 2002). 

Radio emission is generated in the unipolar-inductor 
model by electron cyclotron maser emission from current- 
carrying electrons in the Io flux tube. The free energy for 
maser emission is provided by a loss-cone instability in 
the reflected electron velocity distribution. The observed 
anti-correlation between infrared footpoint emission and Io- 
controlled Jovian decametric radiation supports the hypoth- 
esis that the radio emission is driven by reflected electrons 
(Connerney et al. 1993). In addition, the negative frequency 
drifts evident in the fine-frequency structure of Jovian de- 
cametric radiation (the Jovian S-bursts) is consistent with 
maser emission driven by reflected electrons (Ellis 1974). 
Theoretical models based on electron-cyclotron maser emis- 
sion can explain the observed high brightness temperatures 
(>10 17 K, Dulk 1970), 100% elliptical polarization (Dulk, 
Lecacheux & Leblanc 1992) and the radiation beaming pat- 
tern. Ground-based and dual spacecraft observations have 
confirmed the theoretical prediction that Io-controlled emis- 
sion consists of hollow conical radiation beams, with opening 
angles < 80° and beam widths < 2° , with the source loca- 
tion at the foot of the Io flux tube (Dulk 1967; Hewitt et al. 
1981; Maeda & Carr 1992; Kaiser et al. 2000). 

Recently, two ultra-short-period (P < 10 min) X-ray 
sources RX J1914+24 (Cropper et al. 1998; Ramsay et al. 
2001, 2002b) and RX J0806+15 (Israel et al. 2002; Ramsay, 
Hakala & Cropper 2002a) were discovered. The remarkable 
characteristics of these sources are that there is only a sin- 
gle period for both the X-ray and optical emission, which is 
interpreted as the orbital period (Cropper et al. 1998) and 
the X-ray and optical emission are anti-phased (see Ram- 
say et al. 2001). The inferred orbital periods are too short 
to be from a binary system containing a white dwarf and a 
main-sequence star. Consequently, these X-ray sources are 
proposed to be compact binaries containing two white-dwarf 
stars, which may be in synchronous rotation with the orbit 
(Cropper et al. 1998). Their X-ray and optical/IR variations 
are contrary to those of typical accretion-powered sources, 
which often show multiple periods, corresponding to the or- 
bital rotation, spin of the accreting star and beats. Moreover, 
the absence of optical polarization and helium emission lines 



does not support the scenario involving mass transfer and 
accretion. 

Wu et al. (2002) proposed a unipolar-inductor model 
to explain the peculiar X-ray and optical/IR variations of 
RX J1914+24, based on the unipolar-inductor model for 
the Jupiter-Io system. In the unipolar-inductor model for 
white-dwarf pairs, Jupiter is replaced by a magnetic white 
dwarf, and Io is replaced by a nonmagnetic white dwarf. 
The X-ray emission is powered by resistive heating of the 
magnetic white-dwarf atmosphere at the footpoints of the 
magnetic flux-tubes carrying the electric currents flowing be- 
tween the white dwarfs. This is analogous to the ultraviolet 
footpoint emission of the Galilean moons. The anti-phased 
optical/infrared emission emanates from the irradiated face 
of the nonmagnetic white dwarf, which is visible when the 
X-ray source is directed away from the observer. The or- 
bital period gradually decreases in the unipolar-inductor 
model because the orbital evolution is driven by power losses 
through gravitational radiation, and is therefore consistent 
with timing measurements of the X-ray pulses from RX 
J1914+24, in which the orbital period decreases at a rate 
consistent with gravitational radiation losses (Strohmayer 
2002). This is in contrast to models for the X-ray emission 
involving mass transfer from one white dwarf to another, 
which predict an increase in the orbital period. In contrast 
to the Jupiter-Io system, plasma-inertia effects are unlikely 
to be significant in the white-dwarf context because of the 
relatively high magnetic field strengths and low plasma den- 
sities (Wu et al. 2002). 

It has also been proposed that unipolar induction can 
operate in isolated magnetic white-dwarf systems with ter- 
restrial planets (Li, Ferrario & Wickramasinghe 1998). In 
this case the planet plays the role of the non-magnetic white 
dwarf in magnetic-nonmagnetic white-dwarf binary systems. 
The induced e.m.f. and dissipated power in the white-dwarf 
systems is significantly higher than in the planet-moon sys- 
tems because the magnetic moment is higher for white 
dwarfs than for planets, and both white dwarfs and plan- 
ets are typically larger than planetary moons. 

We predict here, by direct analogy with the Jupiter- 
moon systems, that compact magnetic-nonmagnetic white- 
dwarf pairs and planets orbiting magnetic white dwarfs are 
strong electron-cyclotron maser sources. The large e.m.f. 
produced by the orbital motion of the system supplies 
high-energy electrons, and the (dipolar) geometry of the 
white-dwarf magnetic field and the low-plasma-density en- 
vironment between the two stars or between the star and 
the planet can develop instabilities that generate electron- 
cyclotron maser emission. 

Around 100 extrasolar planets are currently known. 
With the exception of pulsar planets, these planets, which 
have masses comparable to Jupiter and orbit around solar- 
like stars, were detected with the radial-velocity technique. 
Other methods proposed for detecting Jovian (gas-giant) 
planets around white dwarfs include observation of the ex- 
cess thermal infrared emission associated with the Jovian 
planet (Ignace 2001; Burleigh, Clarke & Hodgkin 2002), or 
of re-emitted hydrogen recombination lines from the plan- 
etary atmosphere for close planetary orbits ( < a few AU) 
around a hot white dwarf (Chu et al. 2001). We propose here 
that radio observations can reveal the presence of terrestrial 
planets in close orbits around white dwarf stars. These sys- 
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Figure 1. The unipolar-inductor model for white-dwarf pairs and 
white-dwarf planetary systems, where E is the induced electric 
field across the nonmagnetic white dwarf or planet. 



tems are remnants of main-sequence planetary systems in 
which terrestrial planets survive the stellar expansion phases 
and migrate to a stable close orbit around the white dwarf. 

This paper is organized as follows. In §2 we use a sim- 
ple model to calculate flux densities for maser emission, as- 
suming parameters appropriate in the environment of white- 
dwarf pairs and of white-dwarf planetary systems. A more 
self-consistent consideration which takes into account the 
geometry and properties of the current sheets along the flux 
tubes will be presented in a future paper. In §3 we discuss 
briefly the astrophysical implications of the results obtained 
from our calculations. 



2 MODEL OUTLINE 

In this Section, we outline the model used to predict the 
flux densities for unipolar-inductor-driven radio emission in 
white-dwarf systems. 



2.1 The unipolar-inductor current circuit 

Figure 1 shows a schematic illustration of the unipolar- 
inductor model used in our calculation. For the magnetic- 
nonmagnetic white-dwarf pair, the system parameters are 
the magnetic white dwarf's mass Mi and magnetic moment 
H, the non- magnetic white dwarf's mass M2, and the bi- 
nary period P. The radii of the magnetic and non-magnetic 
white dwarfs, Ri and R2 respectively, are obtained from the 
Hamada & Salpeter (1961) mass-radius relation (assuming 
helium white dwarfs), and the orbital separation a is cal- 
culated from the orbital period and the masses of the two 
white-dwarfs using Kepler's third law. For the white-dwarf 
planetary system, the nonmagnetic white dwarf is replaced 
by a planet, and we use the same notation for the mass 
and radius (i.e. M2 and R2) for the planet. However, the 
planetary radius is treated as a free parameter instead of a 
quantity deduced from a mass-radius relation. 

We define an asynchronism parameter £, which is the 
ratio of the spin angular speed of the magnetic white dwarf 



to the orbital angular speed. In terms of the orbital param- 
eters and the asynchronism parameter, the e.m.f. across the 
nonmagnetic white dwarf is 



-2/3 



(1) 



(see Wu et al. 2002). For white-dwarf pairs with P ~ 10 min 
and Hi ~ 10 32 G cm 3 , a seemingly small departure from 
synchronous rotation (e.g., (1 — £) ~ 0.001) can produce a 
sufficiently high e.m.f. to power X-ray emission at the mag- 
netic footpoints with a luminosity a thousand times the solar 
bolometric luminosity. 

Most of the electrical energy in magnetic-nonmagnetic 
white-dwarf pairs is dissipated in the flux tube footpoints at 
the surface of the magnetic white dwarf. This is because the 
cross-sectional area of the flux tube connecting the two white 
dwarfs at the footpoints is much smaller than the cross- 
sectional area at the nonmagnetic white dwarf. White-dwarf 
planetary systems will also dissipate most energy at the foot- 
points, provided that the conductivity of the planet is higher 
than the white-dwarf atmospheric conductivity. This condi- 
tion is satisfied for a planetary body composed entirely of an 
Earth-like core, after the outer planetary layers have evap- 
orated away at an earlier stage in the system's evolution. 
In the case where a mantle layer is still present, the bulk 
of the electrical power would heat the planet if the mantle 
conductivity is lower than the white-dwarf atmospheric con- 
ductivity (Li et al. 1998). We consider only the case with an 
Earth-like planetary core in this study. 

The effective resistance of the magnetic white dwarf at 
the footpoints of the flux tube is 



Hi 



2(71 i?2 



3/2 



(2) 



(3) 



(see Wu et al. 2002) where ai is the white-dwarf atmospheric 
conductivity, and a is the orbital separation. The current 
flowing in the circuit is 

/=-* 

Hi ' 

and is confined to a current sheet at the surface of the flux 
tube. The number density of current-carrying electrons rifoot 
at the footpoint is given by 

I 

evA ' 



"-foot 



(4) 



where e is the electron charge, v is the average velocity of 
the current-carrying electrons, and A is the area of the arc- 
shaped footpoint. We can approximate A by 



A « TV 3 Rtube(Rl) 



(•») 



where 8 is the relative width of the current sheet, and 
^tube(^i) is the flux-tube radius at the footpoint. Assuming 
a dipolar white-dwarf magnetic field, the flux-tube radius 
varies with distance r from the magnetic white dwarf as 



Rtubc(r) » R 2 (£) 



3/2 



(6) 



It follows that the number density of current-carrying elec- 
trons varies with distance as 

n c (r) ?s rifoot ( — ') ■ (7) 
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Table 1. System parameters: mass M2 and radius R2 of the nonmagnetic white dwarf (first row) and planet 
(second and third rows), period P, asynchronicity parameter f, induced e.m.f. <£, footpoint flux-tube radius 
^tubc(-Ri)i an d number density of current-carrying electrons at footpoint rif oot . 



M 2 (M ) 


i?2 (cm) 


P (hr) 




$ (statvolt) 


-Rtubc(-Ri) (cm) 


"foot (cm 3 ) a 


0.5 


9.6 x 10 8 


0.14 


0.999 


4.0 x 10 4 


2 X 10 7 


1.7 x 10 12 


1.2 x icr 6 


3.5 x 10 s 


10 





960 


1.4 X 10 5 


6.1 x 10 12 


1.2 x icr 6 


3.5 x 10 8 


28 





88 


5.0 x 10 4 


1.6 x 10 12 



assuming 8 = 0.05 rad and v = 0.1c 



2.2 Formation of loss-cone electron distribution 



As the downward-propagating, non-thermal, current- 
carrying electrons reach the polar footpoints, they are lost 
to collisions in the dense white-dwarf atmosphere. The dis- 
sipated electrical power produces the footpoint X-ray emis- 
sion. We assume that collisional losses become significant at 
a height O.Oliii above the white-dwarf surface, correspond- 
ing to the thickness of the white-dwarf atmosphere (e.g. Kip- 
penhahn & Wiegert 1989). 

Assuming adiabatic motion of electrons along the field 
lines, the electron speed is constant, together with the adi- 
abatic invariant sin 2 a/B, where the electron pitch angle 
a is the angle between the velocity v and the magnetic 
field B. Higher up the flux tube (away from the footpoint), 
the precipitating electrons occupy the range of pitch angles 
< a < aic, with the critical angle ai c given by 



ai c = sin 



(B_\ 1/2 ' 



(8) 



where B is the local magnetic field strength and Bf is the 
field strength at the footpoint. Thus the precipitating elec- 
trons map into a decreasing range of pitch angles with in- 
creasing distance from the footpoint. The other downward- 
propagating electrons with pitch angles a > a% c are re- 
flected, before reaching the white-dwarf atmosphere, by a 
magnetic mirror formed by the converging field lines at the 
footpoint. Hence, the distribution of the reflected (upward- 
propagating) electrons at this position has a deficiency in 
electrons with pitch angles a < ai c . This is the loss cone. 
The critical (loss-cone opening) angle Qi c decreases from 
7r/2 rad, at the footpoint, with increasing distance along 
the flux tube. 

The loss-cone distribution is modeled here with the 
Gaussian form: 



f(v,a) = f (v)F(a), 

with the velocity distribution 



fo(v) = 



W3. 



exp 



(27T j ■ f/(,e ri c 

and the angular distribution 



2^ 



f ° 

F(a)={ 1 



[a - (ai c - Aa) 



a < aic — 
a > a ic 
otherwise 



Aa 



(9) 



(10) 



(11) 



where Aa is the loss-cone edge width, ni c is the number 
density of loss-cone electrons, and V\ c = (fcBTic/n'ic) 1 '' 2 is 
the thermal- velocity spread for the loss-cone temperature Ti c 



A V. 




The loss-cone distribution is illustrated as contours 
space, where u 11 and v± are the components of the 



Figure 2 

in — v_ 

velocity v parallel and perpendicular to the magnetic field B. 
The greyscale contours at low velocities correspond to the thermal 
background electron distribution, with number density ra t h and 
temperature T t h <C Ti c . A resonance circle (see Appendix A), 
with centre Do and radius V, is illustrated by a dashed line. The 
angle A<f> corresponds to the segment of the resonance circle which 
intersects the loss-cone edge. 



(fee is Boltzmann's constant). The loss-cone distribution is 
illustrated in Figure 2, as contours in tin — v± space (the com- 
ponents of v parallel and perpendicular to B). We note that 
our calculations assuming an idealized loss-cone distribution 
tend to overestimate the maser growth rate. In practice, the 
maser will operate close to a state of marginal stability due 
to competition between formation and maser-induced relax- 
ation of the loss-cone instability (Robinson 1991). These 
considerations will be included in a more comprehensive 
model in a future paper. 

We adopt the loss-cone instability in our model based on 
the following strong evidence for loss-cone-driven maser pro- 
ducing Io-controlled Jovian decametric radiation: (i) Con- 
nerney et al. (1993) reported an anti-correlation between the 
occurrence of Jovian decametric radiation and the presence 
of an infrared auroral footprint as a function of Jovian lon- 
gitude. Radio emission is preferentially emitted at flux-tube 
footpoint longitudes with higher field strengths, where the 
magnetic mirror is more efficient. The field strength varies 
as a function of footpoint longitude due to the asymme- 
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tries in the tilted and non-centred Jovian magnetic field. 
At these longitudes a higher proportion of electrons are re- 
flected, driving maser emission at the expense of precipi- 
tating electrons, which produce the footpoint emission, (ii) 
ft is well established that the radiation pattern of Jovian 
decametric radiation is a hollow cone, with opening angle 
<80° and beam width < 2° (Dulk 1967; Maeda & Carr 
1992; Kaiser et al. 2000). This is consistent with the theo- 
retical prediction for loss-cone driven masers (Hewitt et al. 
1981). 

We note that an alternative mechanism, where the 
maser instability is driven by a "shell" electron distribu- 
tion, is also proposed to be relevant to astrophysical con- 
texts (Ergun et al. 2000). The shell electron distribution, 
which drives Earth's auroral kilometric radiation, is created 
by magnetic-field-aligned electric fields in the auroral zones. 
We argue that a shell maser is not applicable to Io-controlled 
Jovian radio emission because it produces a thin-disc radi- 
ation pattern, beamed exactly perpendicular to the mag- 
netic field direction. It is unlikely that refraction or scatter- 
ing can transform thin-disc beams into the observed hollow- 
cone beams, which have opening angles 9 < 80° , and angu- 
lar widths < 2° . However, we cannot rule out the opera- 
tion of the shell maser in white-dwarf systems. Nevertheless, 
the higher efficiency of the shell maser will produce higher 
flux densities than predicted by the loss-cone maser model 
adopted in this study, thereby improving the prospects for 
radio observation of white-dwarf systems. 



2.3 Flux density spectra 

Melrose & Dulk (1982) derived semi-quantitative expres- 
sions for the loss-cone-driven electron-cyclotron maser 
growth rate T(lo, 6), emission bandwidth Alo/lo and angular 
spread A8, as summarized in Appendix A. The wave modes 
generated by the electron-cyclotron maser depend most crit- 
ically on the ratio lo p /Q, c in the source region (Hewitt, Mel- 
rose & Ronnmark 1981), where uo v is the electron plasma 
frequency. The escaping wave modes (i.e., modes in which 
waves can escape from the source region to reach a distant 
observer) favoured by the electron-cyclotron maser are the 
fundamental and second-harmonic rc-mode and the funda- 
mental o-mode. The magnetoionic x-mode and o-mode have 
opposite senses of circular polarization. Higher harmonics 
have relatively low growth rates and are unlikely to grow 
to observable levels. In Appendix B we summarize the pa- 
rameter ranges for these three modes in which the electron- 
cyclotron maser operates. 

We estimate the peak brightness temperature in the 
source region by assuming saturation of the reactive version 
of the electron-cyclotron maser instability, i.e., where elec- 
tron phase-bunching effects are significant. This assumption 
is made by analogy to Io-controlled Jovian decametric radio 
emissions, where the highest flux densities are attained by 
the subclass referred to as Jovian S-bursts, which exhibit 
properties consistent with phase-bunching effects (Carr & 
Reyes 1999; Willes 2002). A reactive growth mechanism sat- 
urates at wave amplitudes where electrons become trapped 
in the field of the growing waves. The saturation wave am- 
plitude corresponds to where the trapping frequency lot be- 
comes comparable with the growth rate, such that electrons 



are trapped by the wave before they can contribute to wave 
growth. The trapping frequency is given by 

'eEk\ 1/2 



LOT = 



(12) 



(Melrose 1986, p. 70), where E is the wave electric field, and 
k is the wavenumber, with k = lo/c for waves with refrac- 
tive index n w 1, which is valid for x-mode and o-mode 
waves sufficiently above their respective cutoff frequencies. 
The peak wave electric field i5 m ax can be obtained by equat- 
ing the trapping frequency (Eq. 12) with the maser growth 
rate (Eq. A8), with 



E n 



rn c cF 



It follows that the peak energy density 



F 2 

W = ax 

47T 



m 2 c 2 T 4 
4ire 2 L0 2 



(13) 



(14) 



The relation between brightness temperature Tb and energy 
density W in the source is 



W = k B T h 



(15) 



\2lYcJ \ LO ) 

(Melrose & Dulk 1982), where Alo/lo is the bandwidth and 
AQ is the solid angle filled by the radiation, with Afl — 
2-7TACOS6 1 . Using the derived expressions for the bandwidth 
(Eq. A12) and angular width (Eq. A13) of the loss-cone- 
driven maser, the brightness temperature 

4n 2 c 6 W 
Kbw^o Aa 2 cos 4 a\ c 

where vo is related to the emission angle 6 in Equation (A7). 
The flux density S is related to the brightness temperature 
by 



(17) 



(Dulk 1985), where the integral is taken over the projected 
area of the source for the differential solid angle dQ. As- 
suming a constant brightness temperature over a source size 
equal to the flux tube radius Rtube, this simplifies to 



Ad 2 



(18) 



where d is the distance to the source. 

As the maser emission at harmonic number s, gener- 
ated at a frequency lo = sfl c , propagates away from the 
magnetic white dwarf, it encounters harmonic absorption 
bands at lo « (s + 1) Q c , (s + 2) Q c , . . .. The first absorption 
band is illustrated schematically in Figure 1. We incorpo- 
rate harmonic damping effects in our calculations, using the 
approximate formula for the optical depth of the harmonic 
absorption bands in Appendix C. 



3 RESULTS AND DISCUSSION 

We consider three cases to illustrate the electron-cyclotron 
maser emission in white-dwarf systems. The first is a white- 
dwarf pair, and the other two are white-dwarf planetary sys- 
tems with differing orbital periods. The system parameters 
of these three systems are listed in Table 1. In all cases, the 
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mass of the magnetic white dwarf Mi is fixed to be 0.7 Mq, 
its atmospheric conductivity o\ is 10 14 s _1 , and its magnetic 
moment /i is 10 31 Gem 3 . In the white-dwarf planetary sys- 
tems, the spin of the white dwarf is not synchronized with 
the planetary orbit. We assume here that the spin period 
of the isolated magnetic white dwarf is significantly greater 
than the orbital period, and without loss of generality we set 
£ = 0. The distances to all systems are fixed to be 100 pc. 
For the loss-cone parameters, unless otherwise stated, we 
assume that Aa = 0.05, which corresponds to emission an- 
gular beam widths on the order of a few degrees (by anal- 
ogy with the Jupiter-Io system), the loss-cone temperature 
feTic = 1 keV (where 1 eV = 1.602 x 10~ 12 erg) and the 
footpoint electron- number density n\ c = 10 cm . The loss- 
cone number density varies along the flux tube according to 
Equation (7). For these parameters we are assuming that 
only a very small fraction of current electrons ( < 0.01%) 
contribute to the loss-cone instability. Our flux density esti- 
mates are therefore conservative. 



3.1 Magnetic-nonmagnetic white-dwarf pairs 

3.1.1 Flux density spectra 

Figure 3 shows the peak spectral flux densities for loss- 
cone-driven maser emission from a magnetic-nonmagnetic 
white-dwarf pair over a range of thermal plasma tempera- 
tures T t h and densities nth- The binary system parameters 
are as shown in the first row of Table 1. We will discuss the 
dependence of the peak flux densities on T t h and n t h imme- 
diately below, and subsequently we discuss the dependence 
on the loss-cone parameters (temperature T\ c , electron num- 
ber density n\ c and width Aa). 

At low thermal plasma densities (n th = 10 8 cm -3 ; 
panels a and b), the maser operates in the fundamental 
x-mode over the entire displayed frequency range, up to 
the high-frequency cutoff (above 100 GHz) corresponding 
to the electron cyclotron frequency at the foot of the flux 
tube. Reabsorption of fundamental a;-mode emission at the 
second harmonic layer (see Appendix C) becomes signifi- 
cant as the thermal temperature is increased. Consequently, 
the flux densities at low frequencies are much lower for 
teTth = 10 eV than for fcaTth = 1 eV (compare panels 
b and a in Figure 3). 

At intermediate thermal plasma densities (nth = 10 9 — 
10 cm" ; panels c, d, e and f), fundamental x-mode emis- 
sion is suppressed at low frequencies (see Appendix B). The 
fundamental suppression frequency increases with increas- 
ing thermal plasma density, as indicated by the relative 
positions of the solid arrows in Figure 3, panels e and g. 
Note that there is a frequency overlap between fundamental 
and harmonic x-mode emission because the upper bound 
to harmonic emission is at twice the fundamental x-mode 
suppression frequency. The abrupt drop in flux density at 
high frequencies is due to thermal damping (see Appendix 
A, Equation (A16)). 

At high thermal plasma densities (n t h = 10 11 cm -3 ; 
panels g and h), the fundamental a;-mode emission generated 
above the suppression frequency is strongly absorbed. Only 
the fundamental o-mode and the second harmonic x-mode 
emission generated below (twice) the suppression frequency 



can escape to an observer. No escaping radiation is produced 
below ~ 2.5 GHz, due to suppression of these modes. 

The peak flux densities also vary with the loss-cone tem- 
perature Ti c . Figure 4 shows the flux density spectra for the 
same parameters as Figure 3, except with felic = 10 keV. 
The peak flux densities are consistently lower than in Fig- 
ure 3, typically by four to five orders of magnitude. The de- 
pendence on the loss-cone electron number density is S oc nf c 
(from Eqs. (A8), (14) and (16)). For example, increasing m c 
by a factor of ten, to 10 cm , shifts all the flux density 
curves up by four orders of magnitude in Figures 3 and 4. 
The peak flux densities also decrease with increasing loss- 
cone width Aa. For example, doubling the loss-cone width 
to Aa = 0.1 rad lowers the flux density curves in Figures 3 
and 4 by an order of magnitude. 



3.1.2 Observability 

We discuss here the observability of maser emission from 
white-dwarf pairs by adopting the parameters used in Fig- 
ure 3 as an illustrative example. We consider two observing 
frequencies, 1.465 GHz and 43 GHz, which are available for 
the Very Large Array (VLA) radio telescope. We mark these 
frequencies by grey lines in Figure 3. In a recent search for 
radio emission from Jupiter-like extrasolar planets using the 
VLA, Bastian et al. (2000) considered a sensitivity limit of 
0.02 — 0.07 mjy at 1.465 GHz. We assume a sensitivity limit 
of 0.1 mjy (1 Jy = 10" 23 erg cm' 2 s" 1 Hz" 1 ) at 1.465 GHz 
and 43 GHz in our discussion. This is a conservative estimate 
for the sensitivity at 43 GHz because the VLA sensitivity in- 
creases with observing frequency. 

Maser emission would be detectable at 1.465 GHz for 
the following ranges of parameters. The fundamental x- 
mode emission is at least an order of magnitude above the 
detection threshold in panel a of Figure 3, and hence would 
be observed over the range of parameters: ni c >10 8 cm~ 3 , 
felth ^ 1 eV, and nth < 10 8 cm -3 . At higher thermal densi- 
ties, the fundamental x-mode emission is suppressed at this 
frequency (see panel c), and at higher thermal temperatures, 
the emission is strongly absorbed (see panel b), reducing 
the parameter range for observability to n\ c > 10 10 cm" 3 , 
l<fc B T th <10 eV, and n th <10 8 cm~ 3 . The fundamental 
o-mode emission is detectable in the parameter range: 
nic > 10 10 cm" 3 , fc B T t h < 10 eV, and n th « 10 10 cm" 3 (pan- 
els e and f). 

Note that for a loss-cone number density of n\ c = 
10 10 cm" 3 , corresponding to ~ 1% of current-carrying elec- 
trons in the loss-cone, the flux densities are eight orders of 
magnitude higher than those in Figures 3 and 4 (with n\ c = 
10 8 cm -3 ). This is still a conservative estimate for the frac- 
tion of loss-cone electrons, considering that in the Jupiter- 
Io system the observed anti-correlation between footpoint 
emission and radio emission implies that a more significant 
fraction of the current-carrying electrons contribute to the 
loss-cone instability. No radio emission would be detected at 
1.465 GHz for nth > 10 11 cm" 3 , because all free space modes 
are suppressed at this frequency (panels g and h). 

Maser emission would be detectable at 43 GHz 
over the following ranges of parameters. The fundamen- 
tal x-mode emission is detectable for n\ c > 10 8 cm" 3 , 
fcBTth < 10 eV, n th <10 8 cm" 3 (panels a and b); 
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Figure 3. Peak spectral flux densities for radio emission from magnetic-nonmagnetic white-dwarf pairs (maximized over emission angle), 
with binary system parameters in Table 1. The loss-cone parameters are A;b71 c = 1 keV, n\ c = 10 8 cm -3 , and Aa = 0.05 rad. Two values 
are considered for the thermal plasma temperature, with fcB^th = 1 eV (panels a, c, e and g in the left column) and fceTn, = 10 eV 
(panels b, d, f and h in the right column); four values are considered for the thermal plasma density, with n t h = 10 s cm -3 (panels a 
and b), n t ^ = 10 9 cm -3 (panels c and d), n t h = 10 19 cm -3 (panels e and f), and n t h = 10 11 cm~ 3 (panels g and h). Spectral density 
profiles for the fundamental x-mode (solid line), the fundamental o-mode (dashed line), and the second harmonic x-mode (dotted line) 
are displayed. The vertical grey bars denote VLA observation frequencies at 1.465 GHz and 43 GHz. 



Tiic^lU" cm", fc B 7th<l eV, 10 8 <n th <10 9 cm" 3 (panel 
c); n Xc > 10 9 cm" 3 , 1 < fc B Tth < 10 eV, 10 8 < n th < 10 9 cm" 3 
(panel d); and n\ c > 10 9 cm -3 , ksTth < 1 eV, 
10 9 <n th < 10 10 cm" 3 (panel e). 

Second-harmonic a;-mode emission is only relevant in a 
small region of parameter space (for n t h ~ 10 11 cm -3 in pan- 
els g and h, in the frequency range 13 — 26 GHz. A large frac- 
tion ( > 50%) of current-carrying electrons must contribute 



to the loss cone (corresponding to n\ c > 10 12 cm 



3 ) in order 

for second harmonic emission to be detected. We therefore 
conclude that second harmonic maser emission only plays a 
marginal role in radio emission from white dwarf systems. 

The predicted radio flux at the Earth from white-dwarf 
pairs is variable, with zero flux when the emission cone is not 
directed towards the observer. For emission cone widths of 
several degrees, the duty cycle is on the order of one percent. 
Thus, for the white-dwarf pair parameters assumed above, 
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Figure 4. Peak spectral flux densities for radio emission from magnetic-nonmagnetic white-dwarf pairs for same parameters as Figure 3, 
except with a higher loss-cone plasma temperature: k^T\ c = fO keV. 



with P — 500 s, the radio pulses have durations of > 5 
seconds. The interval between pulses is equal to the orbital 
period P for the extreme viewing geometry where one edge 
of the emission cone grazes the line-of-sight at one point 
in the orbit. More generally, an emission cone with opening 
angle 6 will intersect the line-of-sight at two points per orbit, 
and within each orbital cycle the two pulses are separated by 
an interval with a value between zero and PO/tt, depending 
on the viewing geometry of the system. 



3.1.3 Unipolar inductor lifetimes 

An important constraint to the probability of detection for 
white-dwarf pairs is the lifetime of the unipolar-inductor 
phase and the binary lifetime. The duration of the unipolar- 
inductor phase is related to the rate of synchronization be- 
tween the magnetic white-dwarf spin and the orbital rota- 
tion, due to Lorentz torques associated with the energy dis- 
sipation in the current circuit (see Wu et al. 2002). For the 
parameters assumed here, this timescale is ~ 10 4 yr. The bi- 
nary lifetime is determined by gravitational radiation power 
losses and is ~ 10 6 yr. Hence, the probability of detect- 
ing close white-dwarf pairs in the unipolar-inductor phase is 
~ 1%. Despite the low probability of detection, two systems 
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RX J1914+24 (Cropper et al. 1998; Ramsay et al. 2001, 
2002b) and RX J0806+15 (Israel et al. 2002; Ramsay et al. 
2002a) with ultra-short periods (P < 10 min), are candidate 
unipolar-inductor white-dwarf pairs and are also potential 
maser sources. If these systems are indeed ultra-short pe- 
riod white-dwarf pairs powered by unipolar induction, their 
discoveries suggest that there is either an abundant pop- 
ulation of such close white-dwarf pairs in the Galaxy, or 
that the unipolar-inductor phase may be longer than the 
above estimate of ~ 1% of the system lifetime. The latter 
possibility may be due to a close white-dwarf pair entering 
the unipolar-inductor phase multiple times; for instance, if a 
(post-unipolar-inductor) synchronous system undergoes in- 
termittent mass exchange that breaks the synchronism and 
initiates a new unipolar-inductor phase. 

3.2 White-dwarf planetary systems 

White-dwarf planetary systems consisting of a magnetic 
white dwarf and an Earth-like core with orbital separation 
of 1 AU (~ 2 x 10 4 i?i) have orbital periods ~ 1 yr. The com- 
ponents in these systems are too far apart for the induced 
e.m.f. across the planetary core to be sufficiently large to ac- 
celerate fast (keV) electrons. Masers operate more efficiently 
in more compact systems which are more effective unipolar 
inductors. Here we consider only compact white-dwarf plan- 
etary systems with orbital periods P w 10 — 30 hr (see Ta- 
ble 1). We first discuss the properties of the maser emission 
from these systems and then other aspects, such as the for- 
mation scenario and lifespan of close white-dwarf planetary 
systems. 

3.2.1 Emission properties 

The principle of maser operation in white-dwarf planetary 
systems is very similar to the white-dwarf pairs discussed 
above. The dependence of the emission properties on the 
loss-cone parameters are qualitatively the same, after rescal- 
ing the system parameters (period and masses) accordingly. 
For systems comprising an Earth-sized planetary core orbit- 
ing a non-spinning magnetic white dwarf, an orbital period 
of ~ 10 hr corresponds to an orbital separation of ~ 200 Ri 
(~ 2 Rq). As shown in Table 1, the induced e.m.f. in these 
systems is about 40 times smaller than for the close white- 
dwarf pairs discussed in § 3.1 but is comparable to the in- 
duced e.m.f. in the Jupiter-Io circuit, which accelerates keV 
electrons. Moreover, because the size of the planetary core 
is two orders of magnitude smaller than the size of the non- 
magnetic white dwarf, the footpoint radius of the current 
flux tubes, and hence the maser source size, is greatly re- 
duced in the white-dwarf planetary systems. 

We show in Figure 5 the flux density spectra for maser 
emission from a white-dwarf planetary system with P = 
10 hr. The predicted flux densities are typically four orders 
of magnitude lower than for the white-dwarf pair systems 
in § 3.1. The parameter ranges for which radio emission is 
detectable with the VLA are reduced from those quoted in 
§ 3.1.2, with the minimum values of n\ c multiplied by a factor 
of ten (which increases the peak flux densities by four orders 
of magnitude). 

The peak flux densities decrease as the orbital separa- 



tion increases. For instance, by doubling the orbital separa- 
tion (corresponding to the system parameters in the third 
row of Table 1), the induced e.m.f. is an order of magni- 
tude lower, but is comparable to the induced e.m.f. in the 
Jupiter- Callisto circuit, which generates radio emission of 
comparable power to the Jupiter-Io circuit (Menietti et al. 
2001). The peak flux densities (not shown) are, on average, 
an order of magnitude lower than in Figure 5; again, the 
decrease is predominantly due to the smaller source size at 
the foot of the flux tube. 

3.2.2 Lifetimes and formation scenario 

Whether or not we can detect maser emission from white 
dwarf planetary systems depends not only on the emission 
properties of the systems but also the number of systems and 
the proportion of the planetary-system lifespan occupied by 
the unipolar-inductor phase. The lifetime of the unipolar- 
inductor phase is limited by the inward drift of the planet 
due to Lorentz torques. The inward drift speed is estimated 
to be (Li et al. 1998) 

4$Ja 2 ,.,„, 
ViMt = -GKhM 2 - (19) 

The characteristic lifetime of the white-dwarf planetary sys- 
tem, r ~ a/f drift, scales with binary period as r oc P 5 and 
with magnetic moment as r oc fi~ 2 ■ The 10-hr and 28-hr pe- 
riod systems considered above thus have lifetimes of ~ 10 7 
and 10 9 yr respectively. In general, the probability of observ- 
ing a white dwarf planetary system in the unipolar-inductor 
phase increases with the binary period P. This trend is offset 
by the decrease in the induced e.m.f and maser flux density 
with increasing P. Increasing the white dwarf magnetic mo- 
ment increases the induced e.m.f and footpoint dissipation, 
but simultaneously reduces the unipolar inductor lifetime. 
For example, if we increase the magnetic moment by a fac- 
tor of 10, to fii — 10 32 Gem 3 for the 28 hour system, the 
lifetime will decrease from 10 9 yr to 10 7 yr. 

There is firm observational evidence for the existence 
of double degenerate white-dwarf systems and of isolated 
magnetic white dwarfs. Given that the low-mass compan- 
ion stars in AM Her type binaries and in intermediate po- 
lars (see e.g. Cropper 1990; Warner 1995) will eventually 
become degenerate, the existence of magnetic-nonmagnetic 
white-dwarf pairs is very likely. However, the existence of 
white-dwarf planetary systems is far less certain. We now 
briefly discuss whether solar-like stellar planetary systems 
could evolve into the maser-emitting white-dwarf planetary 
systems investigated theoretically here. 

Solar-type stars undergo two major expansions towards 
the end of their lives — the first occurs as they reach 
the Red Giant Branch (RGB), and the second occurs on 
the Asymptotic Giant Branch (AGB) at the end of the 
core helium burning. The AGB phase is characterized by 
multiple thermal pulses caused by nuclear instabilities in 
the helium-burning shell. In each thermal pulse, the stellar 
radius expands by approximately 25% and contracts over 
short timescales (Sackmann et al. 1993). Following the AGB 
phase, the star rapidly contracts and eventually becomes a 
white dwarf. 

The expansion phases of a solar-type star have crucial 
effects on the fate of their orbiting planets. There are three 
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Figure 5. Peak spectral flux densities for radio emission from an Earth-sized planetary core orbiting a magnetic white dwarf (maximized 
over emission angle) for two values of thermal plasma temperature (fee^th = 1 eV for panels a, c, e, g; fcrj^th = 10 e V for panels b, d, f, 
h), and four values of thermal plasma density (n t jj = 10 s cm -3 for panels a and b; n t h = 10 9 cm -3 for panels c and d; n t h = 10 10 cm -3 
for panels e and f; nth = 10 11 cm~ 3 for panels g and h). The orbital parameters are shown in second row in Table 1. The loss-cone 
parameters are k^T\ c = 1 keV, n\ c = 10 8 cm — 3 and Act = 0.05. Spectral density profiles for fundamental rr-mode (solid line), fundamental 
o-mode (dashed line), and the second harmonic i-mode (dotted line) are displayed. 



likely scenarios: (i) Engulfment: If the expanding stellar ra- 
dius exceeds the orbital separation during the RGB phase 
or the AGB phase (excluding thermal pulses), the planet 
will evaporate, (ii) Temporary engulfment: This is when the 
stellar radius intermittently exceeds the orbital separation 
during the thermal pulses in the AGB phase. The planet 
lies within the stellar envelope for a short period (between 
10 2 and 10 3 yr) and is only partially evaporated (Sackmann 
et al. 1993; Rybicki & Denis 2001). During the temporary 
engulfment, the dominant force is the "bow-shock" drag, 



which causes the planet to spiral inwards, (iii) Otherwise, 
two competing effects determine the evolution of the plane- 
tary orbits. The first is that the orbital separation increases 
due to mass loss from the star in the giant phase (Sack- 
mann et al. 1993). The second effect is that tidal drag forces 
(through which the orbital angular momentum of the planet 
is transferred to the stellar envelope) become significant as 
the stellar radius approaches the orbital separation. The 
tidal interaction acts to decrease the orbital separation, and 
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it dominates over stellar-mass-loss effects when the stellar 
radius expands beyond a threshold value. 

A prediction for the fate of solar system planets, based 
on various existing solar evolution models (Rybicki & De- 
nis 2001), is that the planets at a greater distance from the 
Sun than Mars will all migrate out due to stellar-mass-loss 
effects. Mars is beyond reach of the expanded red-giant en- 
velope of the Sun, and the tidal interaction is negligible; so 
it will not spiral in. Mercury and Venus will spiral in and 
evaporate. The fate of the Earth is less certain, depending 
on the assumed solar evolutionary model. Rybicki & Denis 
(2001) estimated the spiral-in distance for the Earth in the 
temporary engulfment scenario to be (10 — 70) Rq per ther- 
mal pulse, and predicted that the Earth can survive several 
thermal pulses before being permanently engulfed. The sur- 
vival probability increases with the assumed mass-loss rate 
during the RGB phase in the evolutionary model because 
the number of thermal pulses in the AGB phase is reduced 
(Sackmann et al. 1993). Clearly, even with uncertainties in 
the precise evolutionary behaviour of the star, there can still 
exist a range of initial orbital separations from which terres- 
trial planets migrate to very small final orbital separations 
at the end of the AGB phase, through a combination of 
temporary engulfment and tidal interaction effects. 

According to this scenario, the parameter space for the 
formation of close white-dwarf planetary systems is nar- 
row, and we should not anticipate a large fraction of main- 
sequence planetary systems with terrestrial planets evolving 
to the compact white-dwarf planetary systems that allows 
efficient unipolar induction to occur, producing strong maser 
emission. The evolutionary scenario discussed above has not 
taken into account the drift-in effects due to the Lorentz 
torque (see Li et al. 1998). Nevertheless, detailed modeling 
of the planet-star interaction with correct treatments of stel- 
lar evolution, hydrodynamics and magnetic interaction are 
necessary to estimate this fraction. 



4 SUMMARY 

We propose that magnetic-nonmagnetic white dwarf pairs 
with short orbital periods < 10 min, and planetary systems 
comprising a magnetic white dwarf and an Earth-sized plan- 
etary core with orbital periods < 30 hr are strong, periodi- 
cally variable, and highly polarized radio sources. The radio 
emission is generated by electron cyclotron maser emission 
in a unipolar-inductor current circuit driven by the asyn- 
chronicity between the spin of the magnetic white dwarf 
and the orbital motion of the nonmagnetic white dwarf or 
planet. This model is a direct analogy to Jovian radio emis- 
sions generated in conjunction with the Galilean moons. 

We apply a simple model to predict the flux densities of 
the radio emission for the two polarization modes received 
at the Earth from white-dwarf maser systems, assuming a 
source distance of 100 pc. Other model assumptions include 
a dipolar magnetic field for the magnetic white dwarf, and 
a constant thermal plasma density in the immediate vicin- 
ity of the binary system. The free parameters in the model 
are the binary system parameters (component masses and 
radii, orbital distance and period), the thermal plasma pa- 
rameters (density and temperature) , and loss-cone parame- 
ters (density, temperature, and loss-cone width) for the elec- 



trons which are accelerated in the flux tube current circuit 
connecting the two components. 

Our conclusions for magnetized-unmagnetized white 
dwarf pairs are as follows. (1) White dwarf pairs with orbital 
periods < 10 min are strong radio sources over a broad range 
of free parameters in the model. (2) The frequency range of 
radio emission in various wave modes is determined by the 
primary white dwarf's magnetic moment and the thermal 
plasma density. (3) The radio emission spectrum can be dis- 
continuous, and the two polarization modes (with opposite 
senses of circular polarization) may be detected at different 
frequencies. (4) The detection probability is limited by the 
lifetime of the unipolar inductor phase relative to the binary 
lifetime. (5) The predicted radio emission is periodic follow- 
ing the orbital motion of the system, with the duty cycle 
depending on the width of the hollow-cone emission beam. 

Our conclusions for white-dwarf planetary systems are 
as follows. (1) Earth-sized planetary cores orbiting mag- 
netic white dwarfs with orbital periods < 30 hr are strong 
radio sources, but over a more limited range of parame- 
ters than for white-dwarf pairs. (2) Electrical dissipation 
in the unipolar inductor circuit causes the planet to grad- 
ually drift in towards the white dwarf. White-dwarf plan- 
etary systems have a significantly longer unipolar-inductor 
phase than white-dwarf pairs. (3) The longer orbital period 
for white-dwarf planetary systems (in comparison to white 
dwarf pairs) implies longer pulse times: the pulse timescale 
is about tens of minutes, with pulse separation about tens 
of hours (compared with the pulse timescales of a few sec- 
onds and pulse separation of a few minutes for white-dwarf 
pairs) . (4) Only a small fraction of main-sequence-star plan- 
etary systems containing terrestrial planets evolve to close- 
orbit white-dwarf planetary systems which produce maser 
emission. Nevertheless, the large number of galactic white 
dwarfs increases the detection probability. 
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APPENDIX A: ELECTRON-CYCLOTRON 
MASER GROWTH RATE 

In a magnetized plasma, the resonance condition between 
waves and electrons is 
s , 



-fie — k\\V\\ 



0, 



7 



(Al) 



where s is the harmonic number, fi e is the electron cyclotron 
frequency, and fey and «y are the components of the wavevec- 
tor and velocity parallel to B. Assuming the semi-relativistic 



approximation (see Wu & Lee 1979), the Lorentz factor sim- 
plifies to 

and the resonance condition traces out a circle in the v\\ —v± 
space (see Fig. 2), with centre on the vy-axis 

2 



fcyC 
v = — . 

and radius 

/.2 „2 

V = c 



fcjjc _ 2(U - Sfi e 
2 



Sfie 



1/2 



(A3) 



(A4) 



The electron-cyclotron maser growth rate for harmonic 
s and wave mode a (with a = — 1 for rc-mode waves and 
a — 1 for o-mode waves) is given by 



r.,a(k) 

with 

A s , CT (p,k) 



pA,j(p,k) S(ui - sfie/7 - fei|U||) (A5) 



, sfi e d , , d , , 



47r 2 e 2 C 2 /3j 
wn <7 (l + T|)a(wn -)/9w 



(A6) 



X CT sin6> + (cos6>-n/3||)T CT 



n<r/3x sin 9 



Js(z)+J' s (z) 



(Melrose 1980; Melrose & Dulk 1982), where T a , K a , and n CT 
are the axial ratio, longitudinal polarization and refractive 
index for wave mode a, J s {z) is the Bessel function of the 
first kind (order s) and J' 3 (z) is its derivative. The argument 
of the Bessel function z — um CT /3_i_ sin#/fi c ; and the velocity 
component = (v±,v\\)/c. The resonance condition 

appears as the argument in the 5-function in Equation (A6), 
and the growth rate is evaluated as a line integral over the 
resonance circle. The conditions for driving a perpendicular- 
driven maser (with df/dp± > on the resonance circle) 
are therefore satisfied for the loss-cone electron distribution 
modeled by Equation (9). 

Melrose and Dulk (1982) derived a semi-quantitative 
approximation for the maximum loss-cone maser growth 
rate, where for a given loss-cone angle a\ c , and resonance 
circle centre vo, the maximum growth rate corresponds to 
the resonance circle which touches the inside edge of the loss 
cone, with radius V = vo sinai c . For x-mode waves suffi- 
ciently above the cutoff frequency of the mode, the refractive 
index n « 1, and the resonance circle centre is related to the 
wave angle 9, with 



vo = c cos 9 . 

The maximum growth rate is 

2-irmlc 2 A Stt7 f (v cos Qi c ) A0(ai c , Aa) 



vo 



Aa 



(A7) 



(A8) 



where A<j> corresponds to the portion of the resonance circle 
which passes through the loss-cone edge, as illustrated in 
Figure 2, with 



A0(ai c , Aa) = 7T — 2 sin 



_i / sin(ai c - Aa) \ 
y sin ai c J 



(A9) 
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and 



TTu>pm c VQ sin ce\ c cos a\ c a s {\ + vo sin a\ c T a /c) 



4um e (l+T 2 ) 



' vo sm a\ c cos a\ c 



v 



(A10) 



where a s = 4s 2s /[2 2s (s!) 2 ] = 1 for s = 1,2. The maximum 
growth rate occurs at frequency a; max ,s, given by 



= s O 1 + 



Vq cos 2 aic 
"2c 2 



(All) 



The bandwidth is obtained by considering the range of 
radii of resonance circles within the loss cone edge, for con- 
stant emission angle (corresponding to a fixed centre of the 
resonance circle). Then, from Equation (A4), the bandwidth 
is 



AuJ v A 2 . „ 

—p— « —, r Aa cos aic sin 2qi c 



(A12) 



This expression generalizes the result obtained by Hewitt et 
al. (1982), for small ai c , to all loss-cone opening angles. The 
angular width A9 is estimated from the range of vo (centre 
of resonance circle) within the loss-cone edge, with 

A cos 6 = — Aa cot ai c • (A13) 

c 

The effect of damping by thermal electrons also needs 
to be considered. Figure 2 shows a resonance circle which 
intersects the uy-axis, with v\\ = vo — V, at the velocity 
v* where the loss-cone distribution merges into the thermal 
background, with 



r = In ( ^-) , (A14) 

assuming that Vth <Si Vi c and n\ c < nth- Those resonance 
circles with 



wo — V < v* , 



(A15) 



sample the strong damping region associated with the ther- 
mal electrons, and the negative contribution to the inte- 
grated growth rate in the thermal region typically far out- 
weighs the positive contribution in the loss-cone region. 
Hence, from Equations (A3), (A4) and (A14), cyclotron- 
maser radiation is strongly damped for resonance circle cen- 
tres less than the critical value 



vo < 



1 — sin aic V V n\ c J 



(A16) 



In this analysis, we approximate the onset of thermal damp- 
ing by a step function in Equation (A16). 



APPENDIX B: MODE SUPPRESSION 

The electron-cyclotron maser can operate in various plasma 
wave modes; however only waves generated in the magne- 
toionic x-mode and o-mode can escape to a distant observer. 
Waves in other modes are trapped in the source region. For 
maser emission to be generated in the escaping modes, the 
peak frequency u max , s (Equation All) must exceed the rel- 
evant cutoff frequency, which for the x-mode is (Melrose, 
1986, p.171) 



and for the o-mode is 



2x1/2 



(Bl) 



(B2) 



With increasing plasma density or decreasing magnetic 
field strength, the ratio u p /Q, e increases such that cj maXjS lies 
below the cutoff frequency, and successive harmonics are 
suppressed. For lj p /Q c < 0.3, the condition w max ,i > ^ x is 
satisfied. In this case, growth of the fundamental x-mode sat- 
urates the maser, because it has a significantly higher growth 
rate than other modes (Melrose, Dulk & Hewitt 1984). For 
higher values of w p /fl e , fundamental a;- mode emission is sup- 
pressed, allowing wave growth in the fundamental o-mode 
and second harmonic x-mode, which have comparable max- 
imum growth rates. The growth rates of higher harmonic 
x-mode and o-mode waves are too low to produce observ- 
able radiation. Trapped modes may also grow (e.g., z-mode, 
Bernstein modes; Melrose, Dulk & Hewitt 1984; Willes & 
Robinson 1996); however, these can only escape to an ob- 
server after conversion to a free-space mode by linear mode 
conversion or nonlinear wave-wave processes (Melrose 1991; 
Willes & Robinson 1996). We only consider growth of free- 
space modes in this paper. Growth of trapped modes may 
lower the effective growth rates of the fundamental o-mode 
and second-harmonic a;-mode assumed in this paper if the 
maser does not saturate (Melrose, Dulk & Hewitt 1984). 

Assuming saturation of the maser, only fundamental x- 
mode emission is generated when the condition cj max ,i > lo x 
is satisfied. Otherwise, fundamental o-mode emission is gen- 
erated in the parameter range ui < w max ,i < w x , and har- 
monic x-mode emission is generated for o> ma x,i < cu x < 
w ma x,2. Hence there is a wide range of w p /O e over which fun- 
damental o-mode and second harmonic x-mode are simulta- 
neously generated. These ranges vary with emission angle 
9, through Equations (A7) and (All). We note that in a 
more self-consistent treatment of maser emission, where the 
maser operates in a state close to marginal stability (Robin- 
son 1991), the fundamental o-mode and second- harmonic 
x-mode compete more effectively with the fundamental x- 
mode over a wider range of parameters. 



APPENDIX C: HARMONIC DAMPING 

As maser emission, emitted at co « sQ , propagates away 
from the source the local cyclotron frequency decreases with 
decreasing magnetic field strength, and the radiation ap- 
proaches a harmonic absorption layer, where to » (s+ l)f2 e - 
Absorption at higher harmonic layers (s + 2, s + 3, . . .) may 
also occur, but the absorption strength falls off rapidly with 
increasing harmonic number. Assuming a slab geometry, 
with linearly varying magnetic field strength (scale length 
Lb ~ Ri/3), and constant thermal velocity Vth, the optical 
depth at the sth layer for x-mode radiation satisfies (Kuncic 
& Robinson 1993) 



T a ,s(oJ,0) = 



■kL-bu s w p /Khsinfl 
c 2 s s\ u) 2 V c 
(l + T„cos6j 2 +TtV t i/c 2 



1+T 2 



(CI) 



